Fractional Quantum Hall State in Coupled Cavities 
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We propose a scheme to realize the fractional quantum Hall system with atoms confined in a 
two-dimensional array of coupled cavities. Our scheme is based on simple optical manipulation 
of atomic internal states and inter-cavity hopping of virtually excited photons. It is shown that 
as well as the fractional quantum Hall system, any system of hard-core bosons on a lattice in the 
presence of an arbitrary Abelian vector potential can be simulated solely by controlling the phases 
of constantly applied lasers. The scheme, for the first time, exploits the core advantage of coupled 
cavity simulations, namely the individual addressability of the components and also brings the gauge 
potential into such simulations as well as the simple optical creation of particles. 



The achievement of trapping ultracold atomic gases in 
a strongly correlated regime has prompted an interest 
in mimicking various condensed matter systems, thereby 
allowing one to tackle such complex systems in unprece- 
dented ways [i|. A major class of simulable systems, 
distinct from the Hubbard model and spin systems, is 
that in a gauge potential, such as the fractional quan- 
tum Hall (FQH) system. The FQH effect arises when 
a two-dimensional (2D) electron gas is in the presence 
of a strong perpendicular magnetic field at a low tem- 
perature. The Hall resistance of such a system exhibits 
plateaus when the Landau filling factor v takes simple 
rational values [2j. The FQH effect at fundamental fill- 
ing factors v = 1/m for odd integers m (even integers for 
bosons) is accounted for by Laughlin's trial wave function 
(in the symmetric gauge) [3] 
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where Zj = Xj +iyj is the 2D position of the jth electron 
in unit of the magnetic length Ib = \/h/eB with B be- 
ing the magnetic field. The elementary excitation of this 
state is a quasihole (quasiparticle) , which has a fractional 
charge +e/m (— e/m) and obeys the anyonic statistics [4j. 
To simulate such a system in trapped atoms, a major 
challenge is to create an artificial magnetic field as the 
atoms in consideration have no real charge. This is done 
with considerable difficulties, for instance, by rapidly ro- 
tating the harmonic trap [5( , by exploiting electromag- 
netically induced transparency [fjj] , or by modulating the 
optical lattice potential [7|, |8j. Additionally, FQH sys- 
tems are also simulable in Josephson junction arrays [9| . 
Recently, coupled cavity arrays (CCAs) [lfj, [ll|, [l2| 
have emerged as a fascinating alternative for simulating 
quantum many-body phenomena, supported by diverse 
technologies, such as microwave stripline resonators, pho- 
tonic crystal defects, microtoroidal cavity arrays, and so 
forth [13, LLJ, lla|. CCAs have complementary advan- 
tages over optical lattices, such as arbitrary many-body 
geometries and individual addressability [l6l |. Recently, 



theoretical works have shown that the Mott-superfluid 
phase transition of p olaritons [lOf 
berg spin chains [l 21 ] can be realized 



and the Heisen- 
in CCAs. These 
works, however, relied only on globally addressing lasers 
and thus could not highlight the key advantage of CCAs, 
namely, the individual addressability in the sense that 
already they can be done similarly or better in optical 
lattices [17|. Moreover, simulating altogether distinct 
classes of systems such as those of itinerant particles in 
a gauge potential still remains open and this will be es- 
pecially arresting if the particles themselves can be cre- 
ated by a purely optical means. In this paper, we bring 
the Abelian gauge potential into the realm of many-body 
simulations using CCAs. We achieve this by actively ex- 
ploiting the individual addressability, which eventually 
enables great versatility which has not been attainable 
in optical lattices. 

As a concrete example, we introduce a way of simu- 
lating FQH systems in CCAs. To be more specific, we 
consider a FQH system of bosonic particles confined in a 
2D square lattice of spacing a in the presence of a per- 
pendicular and uniform artificial magnetic field B. Non- 
interacting free bosons in a single Bloch band are de- 
scribed by the Hamiltonian 
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where Cj denotes a boson annihilation operator at site 
j and $o = h/e is the magnetic flux quantum. The 
summation runs over nearest neighbor pairs. If we take 
a Landau gauge, this Hamiltonian is written as 
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(3) 
where the positions of lattice sites are represented by 
a(px + qy) with a being the spacing of the lattice. 
Here, a = Ba 2 /$0i the number of magnetic flux quanta 
through a lattice cell, plays a crucial role in character- 
izing the energy spectrum, whose self-similar structure 
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FIG. 1: Schematic representation of a two-dimensional array 
of coupled cavities. Each atom is confined in the intersection 
of two orthogonal cavity modes, which are adjusted to have 
different resonant frequencies. 



FIG. 2: Involved atomic levels and transitions. There are 
two independent Raman transitions mediated by an excited 
level |e) between two ground levels |0) and |1), represented by 
superscripts X and Y, respectively. 



is known as the Hofstadter butterfly [18|. In addition 
to this non-interacting Hamiltonian, we consider a hard- 
core interaction between bosons, which limits the number 
of particles that can occupy one site to a maximum of one. 
In this limit, if we also take a continuum limit a«l, the 
Laughlin state (p} is a very accurate variational ground 
state, where the filling factor v corresponds to the ratio 
of the number of bosons to the number of magnetic flux 
quanta. 

In order to realize the above situation, we consider a 
two-dimensional array of coupled cavities each confining 
a single atom with two ground levels, which will be rep- 
resenting an s = i spin. First, we notify that, aside 
from the additional phases, the Hamiltonian ([3]) in the 
hard-core limit corresponds to that of an s = \ spin 
lattice system with XX interaction, where the creation- 
annihilation operation of the zero and one boson states 
is analogous to the spin flip operation of the spin-down 
and up states. This natural realization of the hard-core 
limit is contrary to the case of optical lattices, wherein it 
is achieved in the limit of strong on-site repulsion [7|, [8| . 
Moreover, as will be seen later, every aspect of the system 
is optically controlled: bosons are created by simple opti- 
cal pulses and the phases in the Hamiltonian are adjusted 
simply by controlling the phases of applied lasers. This 
optical control of the system would greatly simplify the 
experiments, compared to the previous schemes involv- 
ing mechanical modulations of the system. Although in 
this work we mainly consider the FQH systems, another 
great advantage is that unlike the previous schemes for 
optical lattices any Abelian vector potential on a lat- 
tice can be also simulated simply by adjusting the laser 
phases in accordance with the formula ([2]). The creation 
of a quasiexcitation, which is achieved by adiabatically 
inserting a flux quantum through an infinitely thin mag- 
netic solenoid piercing the 2D plane [3|, again reduces 
to the matter of adiabatically changing the laser phases 
accordingly. It can be moved along the lattice cells by 
modulating the laser phases, which would be useful for 
testing the fractional statistics. 

Schemes for realizing the spin exchange Hamiltonian 
in an array of coupled cavities have been established 
in recent papers [12J . In these schemes, the spin ex- 
change is mediated by inter-cavity hopping of virtually 



excited cavity photons. An important difference in the 
present case is that the spin exchanges are associated 
with phase changes depending on their locations and di- 
rections, which obviously cannot be excluded by local 
phase transformations. For this reason, we introduce 
an asymmetry in the 2D geometry of coupled cavities, 
as shown schematically in FIG. [1] where two orthogonal 
cavity modes along the x and y directions have different 
resonant frequencies. Realizing this geometry would be 
viable in several promising models for coupled cavities, 
such as photonic bandgap microcavities [13[ and super- 



conducting microwave cavities |14[. We assume the fre- 
quency difference between the two modes is much larger 
than the atom-cavity coupling rates. In this way, either 
direction of the spin exchange can be accessed individ- 
ually by choosing the laser frequency. We note, how- 
ever, that the above asymmetry is, in fact, not essential 
for our purp ose. For example, an array of microtoroidal 
cavities [151 ] . in which case realizing such a geometry is 
nontrivial, can be also used by involving more lasers in 
the scheme. We discuss this point later. 

We consider a simple atomic level and transition 
scheme as shown in FIG. [2] The atom has two ground 
levels |0) and |1), and an excited level e). There are 
two cavity modes along the x and y directions, whose 
annihilation operators are denoted by a x and a Y , re- 
spectively. The atom interacts with these cavity modes 
with coupling rates g x and g , and with detunings A x 
and A Y , respectively. Two classical fields with (complex) 
Rabi frequencies Q, x e~ 10 and fr e arc applied, re- 



spectively, as in the figure. 
Hamiltonian reads 
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where J x (J Y ) denotes the inter-cavity hopping rate 
of the photon along the x (y) direction, and the sub- 
script (p, q) represents the cavity site. As mentioned 
above, we assume A x — A Y 3> g x ,g Y , and also assume 
A M ^S> g^ ^$> O^, J M . This requires the strong atom- 



cavity coupling in that g^ 3> J M . In this regime, the 
atomic excitation is suppressed, and adiabatic elimina- 
tion leads to an effective Hamiltonian 



H = 



E E [^«(i°m°i)j 



E( \ Jtt V 



h.c. 



(5) 



y„ y t 



/■■'/ 



where ^ = (.<) 2 /A^, u^ = g^W/A^, and ct+ = |1) (0|. 
Here, we have ignored the ac Stark shift induced by 
classical fields, which is negligible in our regime (or it 
may be compensated by other lasers). Again, we assume 
<5'' 3> J M S> o-i^", which can be satisfied, along with the 
above condition, when 



ff P/ A M > jn/ g n > ry/AM. 



(6) 



In this regime, the cavity photon is suppressed, and adi- 
abatic elimination can be applied once more. We extend 
the method in Ref. [19( to keep up to the third order 
terms and take only the subspace with no cavity pho- 
ton. The effective Hamiltonian, in the rotating frame, 
can then be derived as 
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where the parameters are chosen such that t = 
J x (uj x /S x ) 2 = J y (uj y /5 y ) 2 . It is easy to sec that this 
Hamiltonian reduces to the Hamiltonian |(3J) if we adjust 
the phases of the classical fields as 
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The FQH Hamiltonian $5§ is obtained if we adjust these 
phases as 9 X = —pq ■ 2-kol and 9 Y = 0. Note that the 
classical fields for 9 Y can be replaced by one global field. 
In order to check the validity of the adiabatic approxi- 
mation from Hamiltonian ([5]) to ([7]) , we have performed a 
direct numerical diagonalization of Hamiltonian (J2J . We 
take a set of parameters 6^/10 = 10cj m = J M , which cor- 
responds to a case where A^/1000 = g<*/100 = W = J». 
To eliminate the edge effects within a limited computa- 
tional capability, we consider a periodic boundary con- 
dition (i.e., a torus). We consider a 4 x 4 lattice with 
a = 1/4 and two bosons, hence four flux quanta in 
total and the filling factor v = 1/2. In view of the 
fact that the cavity photon is suppressed, we restrict 
our calculation to the subspace wherein the maximum 
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FIG. 3: Energy gap from the ground state to the nine lowest 
excited states for a 4 x 4 lattice in the periodic boundary 
condition in the presence of energy shift — e |1) (1| applied at 
two sites chosen evenly. The filling factor is u — 1/2. 



number of excitations in a cavity is limited to one, i.e., 

<<!<<* + «£>L + (I !> (%,,) < 1- Up to the modifica- 
tion due to the torus geometry and a different gauge [201 ] , 
the ground state should be close to the Laughlin state ([I]) 
with m — 2. From our numerical diagonalization, the 
fidelity between the Laughlin state {^2} and the numeri- 
cal ground state \^c) is found to be Fq = K^I^g)! — 
0.976. We note that when the ideal Hamiltonian $5§ is 
diagonalized for the same 4x4 lattice, the fidelity of 
the ground state is found to be 0.989. The fidelity Fg 
converges to this value as 5^ / J^ and J^/u)^ increase, 
which, however, demands more strong atom-cavity cou- 
pling. Note also that the non-unit fidelity is partly due 
to the finite a, which makes the effect of the lattice ge- 
ometry non-negligible. The ground state fidelity of the 
Hamiltonian (J3j> increases close to one as a decreases [8[ . 

In experiments, the ground state could be prepared 
by the adiabatic transformation, in a similar manner 
as in Ref. |8|. To show this, we consider the above 
4x4 lattice system and deliberately add an energy shift 
— e[(|l) (1|) 00 + (|1) (1|) 2 2]' w hich can be done in experi- 
ments by applying lasers at those sites to induce ac Stark 
shifts. When the energy shift e is sufficiently large, the 
ground state is simply |1) |1) 2 2 w i*h all other atoms 
in state |0). In FIG. ©, we plot the energy gap from 
the ground state to the nine lowest excited states with 
respect to the amount of the energy shift e. The degen- 
eracy of the ground state in the absence of the energy 
shift is due to the ambiguity of the center of mass in 
the torus geometry, which disappears in the plane geom- 
etry [20(1 . Aside from this degeneracy, the excited states 
have finite energy gaps which allow an adiabatic transfor- 
mation. From this figure, it is apparent that the Laughlin 
ground state can be prepared by the following procedure: 

(1) Prepare the atoms in state |1) at sites chosen evenly in 
agreement with the filling factor v, with all other atoms 
prepared in state |0). Initially all lasers are turned off; 

(2) Apply lasers at the chosen sites to induce an ac Stark 
shift — e |1) (1|, with e chosen moderately, e.g., as the de- 
sired value of t. This energy shift does not change the 
atomic state; (3) Gradually increase the Rabi frequencies 
il x and fl Y to reach the desired value of t; (4) Gradually 



decrease the energy shift e to zero. 

The quasiexcitation of the Laughlin state is gener- 
ated when one magnetic flux quantum is adiabatically 
inserted through an infinitely thin solenoid piercing the 
2D plane [3[. In the present system, we can choose the 
position of the quasiexcitation inside a lattice cell. Re- 
calling that the vector potential outside a solenoid is 
given by A s = $o/27rr<^, where r is the distance from the 
solenoid and <p is the azimuthal vector, the effect of the 
solenoid can be easily reflected in the phases of Eq. ((8]). 
Generation of the Laughlin state and the existence of 
the fractionally charged quasiexcitation (in the present 
case, fractionally excited bosons) could be examined by 
directly measuring the individual atoms: for example, by 



measuring the pair correlation functions [2l| . Before the 
measurement, one may turn off all lasers so as to isolate 
the state from further evolution and decoherence. 

Although the atomic excitation is highly suppressed, 
the atomic spontaneous decay is yet a prominent source 
of decoherence. If we denote by 7 the spontaneous decay 
rate of an atom, the effective decay rate of the whole sys- 
tem due to the atomic decay is estimated as iVf,7(fi/A) 2 , 
where Nb denotes the total number of bosons in the sys- 
tem (we omit superscript X or Y for simplicity). On 
the other hand, the energy scale t in the Hamiltonian is 
given by J(il/g) 2 . In view of the condition A 3> g, the 
former is still much smaller than the latter for moderate 
Nb if we assume 7 < J. However, since the excitation 
gap is smaller than t, the attainable system size would 
be restricted in the experimental realization. Although 
the effective decay rate is decreased by increasing A, this 
in turn requires more stronger atom-cavity coupling rate 
g so as to satisfy the condition ©. 

Finally, we stress the point that the asymmetric geom- 
etry introduced in FIG. [1] is not essential. That is, when 
the two orthogonal cavity modes have the same resonant 
frequency, one can also obtain the Hamiltonian in the 
following way: we apply lasers with the same frequency, 
say Wi, in every second row so that they produce the 
spin exchange to the x direction, while applying lasers 
with a different frequency ui2 in the other rows, which 
also produce the spin exchange to the x direction. If we 
choose those frequencies so that \u>\ — wa| ~ <f, they do 
not produce the spin exchange to the y direction. In the 
same manner, we apply lasers with frequencies W3 and W4 
in every second column to produce the spin exchange to 
the y direction. By choosing those four frequencies to be 
sufficiently detuned, we can adjust the associated phases 
independently for each pair of the spin exchange. 
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